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Abstract 

The ion-hose (or fast-ion) instability sets limits on the 
allowable vacuum in a long-pulse, high current 
accelerator. Beam-induced ionization of the background 
gas leads to the formation of an ion channel which 
couples to the transverse motion of the beam. The 
instability is studied analytically and numerically for 
several ion frequency distributions. The effects of beam 
envelope oscillations on the growth of the instability will 
be discussed. The saturated non-linear growth of the 
instability is derived analytically and numerically for two 
different ion frequency distributions. 

1 INTRODUCTION 

With the advent of DARHT-2 and its 2 usee pulse 
concern has surfaced over the ion-hose instability that may 
arise from the beam's interaction with the channel created 
from the background gas via collisional ionization [1]. We 
consider a simple model with a T-dependent neutralization 
(T here is defined as the distance back from the beam head 
divided by (3c). The model we use for the channel is really 
more appropriate for a channel that has been preformed 
however a more correct treatment (which is also more 
complicated) produces the same asymptotic result so we 
adopt this one for simplicity. Let y represent the centroid 
position of the beam and \|/ represents the centroid position 
of the ion channel. We consider the case of a smooth 
external focusing force (the case of solenoidal focusing is 
well represented in the asymptotic solutions by replacing 
k the betatron wavenumber by k 12, one half of the 
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cyclotron wavenumber). 

Since we will also follow the non-linear development 
of the instability we choose as a starting point the 
equations used by Buchanan to describe the coupling of a 
beam and channel, each assumed to have a Gaussian spatial 
profile [2]. The beam has Gaussian radius a while the 
channel has Gaussian radius b such that the parameter R Q 
is given by 

R 2 =a 2 +b 2 . [1] 

This parameter results from integrating the force due to the 
beam over the distribution of the channel and vice versa. 
The model is 
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with y = y I R a and \j/ = \]/ / R a . Here £=kpz, x=(0 o t 
and 

k 2 

[5] 
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where x G is the neutralization time of the background gas 
and is approximately given by 

T o = 10- 9 /P,_sec M 

k 2 is the coupling strength given by (I is the beam current 
and I G is ~ 17 kA) 
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and co is the ion (angular) "sloshing" frequency in the 
field of the beam 
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Here M is the ion mass, q the ion charge and c the speed 
of light. 

2 LINEARIZED EQUATIONS 

If both y and \|/ are small compared to R D equations [2] 
and [3] may be linearized as 
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+ y + £x(y -y/) = 



d 2 y 

dx 2 
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We will solve these equations for a "tickler" excitation, 
that is 

dy(0,x) 
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y(0, X )= ¥ (0, X ) = ^^l = 0. [12] 
ox 



where 



dx 



t + a 2 (y,-y) = 



These equations are appropriate for a beam and channel 
[3] system that are characterized by a single betatron and ion 
"slosh" frequency. We are treating the case of solenoidal 
focusing that we assume is dominant compared to the 
focusing provided by the ion channel. Under this condition 



it is a good approximation to neglect the spread in betatron 
frequency that will result due to the non-linearity of the 
beam-channel force which arises from the non-uniform 
spatial profile of the channel. 

However, it is not a good approximation to ignore the 
spread in ion resonance frequencies which arises from the 
non-uniform spatial profile of the beam. To account for 
this spread we use the "spread mass" model [3], We 
modify the model by splitting the channel centroid into 
"filaments" labeled by a subscript X which characterizes the 
frequency of a particular filament. Equation [10] is thus 
modified as 

^r+*M-y) = o. [13] 

The position of the channel centroid is then found by 
averaging the individual positions of the filaments over a 
distribution function 



W = \f{k)V x dX. 
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For numerical work we will use the conventional definition 
and take a "top hat" distribution where 



/(A) = - for 1-0<A <1. 
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For analytic work we will use a Lorentzian distribution 
(and equation [13] with X 2 instead of X): 
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where 8 is the half-width of the distribution and the range 
of X is from — oo to +°°. 

3 EFFECTS OF ENVELOPE 
OSCILLATIONS 

We now wish to investigate the effects of an envelope 
mismatch in the accelerator on the growth of the 
instability. Since the dominant focusing for the beam is 
provided by the solenoidal field, an envelope mismatch 
will result in a beam radius that varies as 



r h =a(l + jUsin2g) 



[17] 



where we have assumed a particular choice of phase for the 
envelope oscillations without loss of generality. Because 
the channel is formed by the beam we can expect that there 
will be a similar variation for the channel radius. Thus the 
ion resonance frequency will be periodically varying in z. 
This is analogous to the case of "stagger tuning" the 
resonant frequency of cavities to detune the beam breakup 
instability. 

We will investigate this effect by averaging over the 
fast C, and x oscillations of both the channel and beam 
centroid positions [4]. 



First we write the factor g(Q as 

g(g) = l + insm2g. 

Equations [9] and [13] then become 
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and 
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Then we may use the Laplace transform method 
(transforming in x to s and back again) along with 
equations [14] and [16] to obtain 
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We now write y(£,x) as 

y(C,x) = A(c,*y ( "~ 
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where A is regarded as a slowly varying amplitude such 
that 
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Treating (0. as a small parameter, averaging equations [19] 
and [21] in £ over 2%, we find after considerable algebra 
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By Laplace transforming equation [24] in x to s and using 
the method of steepest descents we find the asymptotic 
growth rate 

yoc e T{k °' s ° ] [25] 
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The exponential growth given by equations [25] and [26] is 
shown in Fig. 1. 
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Fig. 1. Asymptotic growth as a function of jj. for 8 =.05. 

We see that a small envelope mismatch can significantly 
reduce the linear growth, particularly for a large growth 
rate. 

4 NON-LINEAR DEVELOPMENT 



It is clear from equations [2] through [4] that when the 
beam and channel displacements become of order R the 
beam-channel force falls off significantly as compared to 
the linear approximation used in equations [9] and [10]. 

We now extend equation [13] into the non-linear range 
(for a Lorentzian distribution) as 
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By taking y = A(£, x)e'^~ x ^and = B x{C^ x ) e 
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with A and B both slowly varying we may average 
equations [2] and [27] to obtain (assuming B»A) 
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[28] 
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Equation [29] can be solved iteratively and integrated with 
equation [16] to find B. This result can be used to 
manipulate [28] into the form 



2y/ 



= 



with 



\B\ 2 =S5^-/ 



1 + 



¥ 



¥s 

2 \ 



¥ s 2 j 



[30] 



[31] 



where y/ = \a\ 2 , y = 32<5 3 , and £ = efy. 1 48 , 

s 

the number of e-folds of linear growth. 
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Fig. 2. Channel and beam position vs. x at the end of the 
accelerator from numerical solution. The blue curve is the 
channel centroid while the purple curve is the beam 
centroid. The top hat distribution (equation [15]) was used 
for 0=0.59. 
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Fig. 3. Channel and beam position vs. \ from equations 
[30] and [31] with 8=0.0939. This maximum value of the 
abscissa corresponds to 10 e-folds of linear growth of the 
beam centroid position, as is the case for figure 2. 

5 CONCLUSIONS 

We have shown that envelope oscillations that lead to 
a periodic detuning of the ion resonant frequency 
significantly reduce the linear growth rate of the 
instability. In addition, when the amplitude of the ion 
channel motion becomes of the order of R , the beam- 
channel force falls off significantly from the linear 
approximation. The betatron motion of the beam/channel 
causes a periodic modulation of the ion resonant frequency 
which increases the effective damping of the oscillations. 
This effect leads to the saturation of the beam centroid 
displacement at an amplitude that is of the order of 28\|/ or 
about 10 - 30% of the channel amplitude. 
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